Microscopic Theory of Spin Transfer Torques in Ferromagnets 
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Spin transfer torque (STT), involving both adiabatic and non-adiabatic contributions, provides 
key mechanisms to current-induced phenomena in ferromagnets. While the adiabatic STT has been 
well studied, current understanding of the non-adiabatic STT is based on macroscopic models. By 
computing microscopically the response of conduction electron spins to a time varying and spatially 
inhomogeneous magnetic background, we derive the adiabatic and non-adiabatic STT in a unified 
fashion. Our result justifies the macroscopic theory (Phys. Rev. Lett. 93, 127204 (2004)) and clarifies 
the physical origin of every coefficient. We also provide a benchmark on the validity of the theory 
and explain why a recent experiment shows deviation. 
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Interplay between current and magnetization is the 
central topic of spintronics 111]. When a current flows 
through a ferromagnetic metal, it becomes spin-polarized 
due to local exchange coupling between conduction elec- 
tron spins and magnetic moments. In turn, spin angu- 
lar momentum is transferred from current to magnetiza- 
tion regarding the constraint of spin conservation, which 
is now known as the spin transfer torque (STT) @, Q. 
STT provides key mechanisms to numerous intriguing 
phenomena in ferromagnets, for example, current-driven 
domain wall motion [j, 0] , spin wave excitations 0, 0] , 
etc. In both fundamental physics and application ori- 
ented studies, STT has aroused enormous attention in 
the last two decades d, However, while magnetiza- 
tion dynamics driven by STT is becoming the core issue 
of spintronics, fundamental physics underlying STT is far 
from clear. 

At present, STT is believed to be divided into adi- 
abatic (reactive) and non-adiabatic (dissipative) contri- 
butions. While the former has been derived microscop- 
ically via different approaches EH, the latter has 
only been justified macroscopically through spin conser- 
vation and Galilean invariance SEMI, whose micro- 
scopic origin is under intense debates. In many recent ef- 
forts, microscopic theories have been developed in generic 
ways HI- lU and in specific contexts 1^- 21 1 , but their co- 
efficients do not reach a consensus, meanwhile, some oth- 
ers even cast doubt on the existence of the non-adiabatic 
STT 22]. On the other hand, experimental measure- 
ments of this torque also exhibits various results that are 
not in agreement 23l425||. 

In this Letter, a microscopic derivation of the mag- 
netization dynamics induced by STT is provided. Wc 
first calculate the response of a conduction electron spin 
to a time varying and spatially inhomogeneous magnetic 
background M (r, t), then obtain the non-equilibrium lo- 
cal spin accumulation 8m (perpendicular to M(r,t)) by 
integration over the Brillouin Zone. Due to the exchange 
coupling between itinerant electron spins and background 
magnetization, the back action of current exerting on 



M (r, t) is proportional to 8m x M, where the adiabatic 
and non-adiabatic STT naturally appear on an equal 
footing. Our result justifies the macroscopic model pro- 
posed in Ref. 3] with all coefficients matched exactly. 
At the end of this Letter, a benchmark on the validity 
of our theory is given and we explain why a recent mea- 
surement on narrow domain walls [25j shows deviations 
from theoretical prediction. 

To avoid confusion of symbols, we adopt the conven- 
tion used in Ref. Q. The exchange interaction between 
itinerant electrons and local magnetic moments is gov- 
erned by the Hamiltonian, 



Hex = ^s-M(r,t) 

M. y ' 



(1) 



where s is the (dimensionless) spin of a conduction elec- 
tron, \M(r,t)\ = M s is the saturation magnetization, 
and 5* denotes the magnitude of background spins. The 
coupling strength J ex can be as large as eV in typical 
transition metals, so that if M(r, t) varies slowly in space 
and time, conduction electron spins will follow the back- 
ground profile when the system is in thermal equilibrium, 
which is known as the adiabatic limit. However, when 
an external current is applied to the system, a small 
non-equilibrium spin accumulation 8m transverse to lo- 
cal M(r,t) is induced. It is this 8m that exerts STT to 
the background magnetization. 

To compute 8m, we first study the spin response 
of an individual conduction electron to the background 
M(r,t) when current is applied. From Eq. (JTJ) , we know 
that local spin-up (majority) and spin-down (minority) 
bands are separated by a large gap A = SJ ex = £± — £f, 
and the associated spin wave functions are denoted by 
| t (r,t)) and | 4(r,t)). The electron is described by a 
coherent wave packet centered at (r c , k c ) 2(| 27 1, 



\W) =J d 6 kw(k)e^ r \k}[c a \r(r c ,t)} + c b \i(r c ,t))} (2) 

where w(k) satisfying J dkk\w(k)\ 2 = k c is the profile 
function, \k) is the periodic part of local Bloch function, 
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FIG. 1: (Color online) Eigenstates of Hamiltonian Eq. ([T]) 
form a set of local spin basis, and defines a local frame that 
moves with n = M /Ms. Components of the conduction elec- 
tron spin s (red) in the local frame are denoted by si, S2, 
and S3. In the tangential plane with normal n, we make a 
coordinate transformation from eg and to n and n x n so 
that everything expressed in the new basis is physical. 



and c a , Cb are superposition coefficients. Now |t(r c ,£)) 
and I 4- (r c ,i)) form a set of local spin basis with quan- 
tization axis being n(r c ,t) = M(r c ,i)/M s , therefore we 
can construct a local frame moving with M(r c ,t) where 
the coordinates are labeled by n, eg, and in Fig. [T] 
The electron spin expressed in this local frame reads, 

s = {sx,s 2 ,s 3 } = 77*077 
= {2Rc(c a c*), -2Im(c a c b *), |c a | 2 - \c b \ 2 }, (3) 

where er is a vector of Pauli matrices, and 77 = [c a ,Q,] T 
is regarded as spin wave function in the local basis. 

The equations of motion are obtained from the univer- 
sal Lagrangian L = (W\ihdt — H\W) through variational 
principle [26| . which involves not only dynamics of r c 
and k c , but dynamics between the two (well separated) 
spin bands as well. The latter represents spin evolution 
with respect to the local magnetization M(r,t) and ex- 
hibits fast rotating character due to the large gap A, it 
should be distinguished with the adiabatic dynamics be- 
tween degenerate bands [27|. Because of space-time de- 
pendence of the local spin basis, Berry gauge connections 
are induced in the effective Lagrangian [28[, 
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where £0 = h{£^{k c ) + £±(k c )) denotes the average of 
local band energy; the gap A couples only S3, which re- 
sembles a local Zeeman energy. The Berry connections 
have both spatial component 



A = ih 
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as a vector potential (V 
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and temporal component 



$ = ih 



t\dt\t) 

mm 



mm 



(6) 



as a scalar potential. The local spin wave functions are 
taken to be | f) = [e~^cos|, e 4 ^ sin 



1 2 sin ■ 



e l 2 cos|] T , where 6 



f] T and I |) = 

(r c ,t) and <f> = 
<p(r c ,t) are spherical angles specifying the direction of 
M(r c , t), whose total time derivatives are 9 = r c -V8+dtO 
and <p = r c ■ V0 + dt<f>. Then the Berry connection terms 
can be written as a unified 2x2 matrix, 



r r - A 



$ = - 
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cos tfcp 
- sin 8d> + i6 



(7) 



It is worth mentioning that freedom exists in the choice of 
local spin wave functions, which leads to the gauge free- 
dom of the Berry gauge connections. More graphically, 
a specified set of spin wave functions corresponds to a 
particular choice of local frame in Fig. [TJ thus the rela- 
tive orientation of the local frame can be rotated about 
n by gauge transformations, and is not physical. But 
everything will be expressed in terms of gauge invariant 
quantities at the end. 

Regarding Eq. ([3]), the spin dynamics is obtained 
through SL/Srj = after some manipulations 28], 
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where r ex = A/H is defined as the exchange time. Eq. © 
describes coherent spin dynamics in the local frame mov- 
ing with M(r c ,t). However, spin relaxation as a non- 
coherent process should also be taken into account. In 
real materials spin relaxation is very case dependent, but 
regardless of the underlying mechanism, it adds a term 
— t~( s — s eq) to Eq. ([8]). where r s f is the mean spin- flip 
time and s eq = {0, 0, 1(— 1)} is the local equilibrium spin 
configuration for the majority (minority) band £f(£|). 
Eq. flSJ should be solved numerically in general, but ap- 
proximation can be made upon the following considera- 
tions: the large gap A results in an extremely small r ea . 
(typically of the order of 10~ 14 ~ 10~ 15 s), thus within 
several t cx s the change of magnetization seen by the elec- 
tron is small, i.e., dtM and (r c ■ V)7W is small com- 
pared to l/r ex . To this end, we define two small pa- 
rameters Ei = T ex sm6<j) and £2 = i~ ex 6 which satisfy 
yje\ + e\ = h\M\/(M s A) < 1. On the same time scale, 
the variations of £1 and £2 are even higher order small 
quantities, thus it is a good approximation to treat £1 
and £2 as constants, by which Eq. ((SJ becomes a set of 
first order differential equations with constant coefficient 
matrix. As a result, it can be solved analytically. As 
explained in the supplementary material [28j, given the 
initial condition s = s eq , the solution of Eq. ([5} for the 
majority band reads, 
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si(t) = - " - [£i(cosf + gsint) +e 2 (sint-^cost)] 

1 + C 1 + ? 

S2 (*) = — l j_ f :2 TTTz t ei ( sin t-£ cost )- £2 (cos t + £ sin t)] 



s 3 (t) = 1 



2 ( e T + ^2) [ cos * + £sint] 



(9a) 
(9b) 
(9c) 



r 



where t = t/r^ is the scaled time and £ = T ex /T s f. 
All the time-dependent terms in Eq. © exhibit oscilla- 
tory pattern along with the exponentially decaying factor 
e - ^*. Since r ex is extremely small, those are fast rotating 
terms (also rapidly vanishing controlled by £), compar- 
ing with the slow dynamics of M . Therefore in deriving 
the dynamics of M those time dependencies are averaged 
out, and we have 



<*i> 

(S3) 
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1, 



(10a) 

(10b) 
(10c) 



and the spin of the electron is s = s eq + 5s where 5s = 
(s\)ee + (s2)e l f > . For the minority band, Eq. (fT0|) only 
differs by an overall minus sign. To express 5s in terms of 
gauge invariant quantities, we need to make a coordinate 
transformation which corresponds to a rotation of basis 
in the tangential plane depicted in Fig. [TJ 
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where Q, = \n\. Then we obtain, 
[n x n + £n] 



<9n ^ dn 



1 + £ 2L "' 'at s at 

+ n x (r c -V)n + £(r c -V)n], (12) 

where n = a 4 n+ (r c -V)ra has been used and r c — — 

is the center of mass velocity. The local non-equilibrium 

spin accumulation is obtained by integration 

6m = HB J d£[9 t (£)g t (£)5 Si + %(£) 9i (£)6 Si ] (13) 

where /is is the Bohr magneton, Sf\,\.{£) is the density of 
states, and g-\^{£) represents the distribution function. 
In weak electric field E and zero temperature, we have 
m,i{£) = /o t 4(£) +^o t 4^-fe%^ where / ot4 (£) 
is the Fermi distribution function without electric field 
and To-f-,4, is the relaxation time. Attention should be 



paid that when the mean spin-flip time T s f is assumed to 
be independent of energy, it is equivalent to introduce it 
either in solving the Boltzmann equation or in Eq. ([8]l. 
and we have chosen the latter. Our target now is to relate 
5m to the charge current 



se 



2> t (£)g t (£) 



d£ t 
dk c 



d£i 

dk c 



After some manipulations, the electric field E and to-cj. 
are eliminated, and we arrive at our central result, 
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where P = (n-j- — n^)/(n^ + nj.) is the spin polarization 
with rw n being the electron density of the two bands at 
Fermi surface, and no — [1b J d£ &{£ )(fo^(£) — foi(£)) 
is the local equilibrium spin density whose direction is 
parallel to the magnetization. Eq. (fl"4"|) is exactly the 
same as Eq. (8) in Ref. @, but the above derivation is 
purely microscopic, and the four terms of Eq. (1141) can 
be traced back to the four terms in Eq. (|12l) . 

From Eq. (J), the STT is T = (l/T ex M s )5m x M, 
which should be added to the Landau-Lifshitz-Gilbert 
equation: dM/dt = jH eff x M+(a/M s )M x 8M/dt+ 
T, where 7 is the gyromagnetic ratio, H e , , is the effective 
magnetic field, and a is the Gilbert damping parameter. 
The final form of magnetization dynamics becomes, 
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(15) 



where u = Pj e /j,B /eM s (l + £ 2 ) is the effective electron 
velocity, and r\ = (n /M s )/(l + £ 2 ) is a dimensionless 
factor. The renormalized gyromagnetic ratio and Gilbert 
damping parameter are 
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where the renormalization originates from the first two 
terms of Eq. p^|) (or Eq. (|14p ). they are determined by 
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the local equilibrium spin density no which exists even 
in the absence of current. Eqs. (|T5|) and ([To]) justify the 
results of previous macroscopic theory Q 

Our microscopic derivation rests on two assumptions: 
local equilibrium can be defined; and M is nearly con- 
stant within several r ex s. The former requires diffusive 
transport which is usually the case in permalloy mate- 
rials; the latter, however, is only true when the charac- 
teristic length of the texture I (e.g., domain wall width) 
satisfies I 3> Vp T ex where vp is the fermi velocity, other- 
wise the solution Eqs. © and (ITU1) are invalid. In a re- 
cent experiment 25| , people measured the non-adiabatic 
torque on very narrow domain walls (1 ~ lOnm) and 
find disagreement with Eq. (|15j) . A rough estimate using 
Vp ~ 3 x 10 5 m/s and A ~ leV tells us that vpT ex is of 
the order of a few angstroms, thus a domain wall of a few 
nm wide cannot be considered I 3> vpT ex . In that case, 
our local solution is no longer a good approximation, and 
STT may exhibits non-local behavior and also oscillatory 
patterns as the time-dependent terms in Eq. ([9]) become 
important. Two final remarks: spin-orbital interactions 
may contribute an anomalous velocity to r c , which brings 
additional terms into Eq. (|14[) ; the spin motive force [291 ] 



E S MF = 7T n ' ( d t n X Vn )> 

2e 

though small, should be considered in a strict sense, and 
the electric field should be replaced by the effective field 
E e ff = E + E S mf in deriving Eq. (TTJ} from Eqs. (|T2|) 
and (fT3|) . This creates an additional contribution to the 
Gilbert damping parameter, which has been studied re- 
cently via quite a different route [3oT| . 

In summary, we have derived the current-induced STT 
including both adiabatic and non-adiabatic contributions 
through a microscopic point of view, which justifies the 
previous macroscopic theory in Ref. 3] with all coeffi- 
cients matched exactly. The validity of our theory is chal- 
lenged when the characteristic length of the background 
magnetization becomes comparable to vjrr ex , which has 
been observed on narrow domain walls [25j • 
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